Ferromagnetism and Metal-Insulator Transition in the Disordered Hubbard Model 
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A detailed study of the paramagnetic to ferromagnetic phase transition in the one-band Hubbard 
model in the presence of binary alloy disorder is presented. The influence of the disorder (with 
concentration x and 1 — x of the two alloy ions) on the Curie temperature T c is found to depend 
strongly on electron density n. While at high densities, n > x, the disorder always reduces T c , at 
low densities, n < x, the disorder can even enhance T c if the interaction is strong enough. At the 
particular density n = x (i. e. not necessarily at half filling) the interplay between disorder-induced 
band splitting and correlation induced Mott transition gives rise to a new type of metal-insulator 
transition. 
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In correlated electron materials it is a rule rather than 
an exception that the electrons, apart from strong inter- 
actions, are also subject to disorder. The disorder may 
result from non-stoichiometric composition, as obtained, 
for example, by doping of manganites (Lai_ x Sr x Mn03) 
and cuprates (Lai^Sr-rCuC^) jjjj, or in the disulfides 
Coi_2;Fe x S2 and Nii-^Coa^ §7 In the first two ex- 
amples, the Sr ions create different potentials in their 
vicinity which affect the correlated d electrons/holes. In 
the second set of examples, two different transition metal 
ions are located at random positions, creating two differ- 
ent atomic levels for the correlated d electrons. In both 
cases the random positions of different ions break the 
translational invariance of the lattice, and the number of 
d electrons/holes varies. As the composition changes so 
does the randomness, with x = or x = 1 corresponding 
to the pure cases. With changing composition the sys- 
tem can undergo various phase transitions. For example, 
FeS2 is a pure band insulator which becomes a disordered 
metal when alloyed with C0S2, resulting in Coi_ K Fe x S2. 
This system has a ferromagnetic ground state for a wide 
range of x with a maximal Curie temperature T c of 120 
K. On the other hand, when C0S2 (a metallic ferromag- 
net) is alloyed with MS2 to make Nii_ x Co x S2, the Curie 
temperature is suppressed and the end compound NiSa 
is a Mott-Hubbard antiferromagnetic insulator with Neel 
temperature Tjy — 40 K. 

Our theoretical understanding of systems with strong 
interactions and disorder is far from complete. For exam- 
ple, it was realized only recently that in gapless fermionic 
systems the soft modes couple to order parameter fluctu- 
ations, leading to different critical behavior in the pure 
and the disordered cases |3) . A powerful method for the- 
oretical studies of strongly correlated electron systems 
is the dynamical mean- field theory (DMFT) |, § §. 
The DMFT is a comprehensive, conserving, and ther- 
modynamically consistent approximation scheme which 



emerged from the infinite dimensional limit of fermionic 
lattice models JtJ. During the last ten years the DMFT 
has been extensively employed to study the properties of 
correlated electronic lattice models. Recently the com- 
bination of DMFT with conventional electron structure 
theory in the local density approximation (LDA) has pro- 
vided a novel computational tool, LDA+DMFT [§, §, 
for the realistic investigation of materials with strongly 
correlated electrons, e.g. itinerant ferromagnets JTo[ . 

The interplay between local disorder and electronic 
correlations can also be investigated within DMFT 
0, H HI, U, ||. Although effects due to coherent 
backscattering cannot be studied in this way [ fTT| , since 
the disorder is treated on the level of the coherent poten- 
tial approximation Jl6[ |, there are still important phys- 
ical effects remaining. In particular, electron localiza- 
tion, and a disorder-induced metal-insulator transition 
(MIT), can be caused by alloy-band splitting. In the 
present paper we study the influence of disorder on the 
ferromagnetic phase. We will show that in a correlated 
system with binary-alloy disorder the Curie temperature 
depends non-trivially on the band filling. In the disor- 
dered one-band Hubbard model we find that for a certain 
band filling (density) n = N e /N a , where N e (N a ) is the 
number of electrons (lattice sites), disorder can weakly 
increase the Curie temperature provided the interaction 
is strong enough. A simple physical argument for this 
behavior is presented. We also find that at special band 
fillings n 7^ 1 the system can undergo a new type of 
Mott-Hubbard MIT upon increase of disorder and/or in- 
teraction. 

In the following we will study itinerant electron fer- 
romagnetism in disordered systems, modeled by the 
Anderson-Hubbard Hamiltonian with on-site disorder 
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where t{j is the hopping matrix element and U is the 
local Coulomb interaction. The disorder is represented 
by the ionic energies e,, which are random variables. 
We consider binary alloy disorder where the ionic en- 
ergy is distributed according to the probability density 
F(e) = x5(e + A/2) + (1 - x)5{e - A/2). Here A is the 
energy difference between the two ionic energies, provid- 
ing a measure of the disorder strength, while x and 1 — x 
are the concentrations of the two alloy ions. For A ^> B, 
where B is the band-width, it is known that binary al- 
loy disorder causes a band splitting in every dimension 
d > 1, with the number of states in each alloy subband 
equal to 2xN a and 2(1 — x)N a , respectively jl6| . 

We solve (0) within DMFT. The local nature of the 
theory implies that short-range order in position space is 
missing. However, all dynamical correlations due to the 
local interaction are fully taken into account. 

In the DMFT scheme the local Green function G an 
is given by the bare density of states (DOS) N°(e) and 
the local self-energy T,„ n as G an — J deN° (e) / (iuj n + — 
Sen — e). Here the subscript n refers to the Matsubara 
frequency icu n = i(2n + 1)tt//3 for the temperature T, 
with [3 = 1/ksT, and /i is the chemical potential. Within 
DMFT the local Green function G an is determined self- 
consistently by 
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together with the k-integrated Dyson equation C/~„ = 
G~ 7 \ + £ CTn . The single-site action Ai for a site with the 
ionic energy = ±A/2 has the form 

U f 

~Tz2 dr <( T )c<r(T-)cl (T (r)c_ (T (T),(3) 
z a Jo 

where we used a mixed time/frequency convention for 
Grassmann variables c CT , c*. Averages over the disorder 
are obtained by (• • -) dis = J deP(e)(- ■ •). 

Since an asymmetric DOS is known to stabilize fer- 
romagnetism in the one-band Hubbard model for mod- 
erate values of U @ || we use the DOS of the 
fcc-lattice in infinite dimensions, N°(e) — exp[— (1 + 

\/2e)/2]/y / 7r(l + \/2e) @. This DOS has a square root 

singularity at e = —l/s/2 and vanishes exponentially for 
e — ► oo. In the following the second moment of the DOS, 
W, is used as the energy scale and is normalized to unity 
PH . The one-particle Green function in Eq. (||) is deter- 
mined by solving the DMFT equations iteratively p"7| , pl| 
using Quantum Monte-Carlo (QMC) simulations p2[ . 
Curie temperatures are obtained by the divergence of the 
homogeneous magnetic susceptibility [pTl p3|. 
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FIG. 1: Curie temperature T c as a function of disorder 
strength A for band filling n larger [panel (a)] and smaller 
[panel (b)] than the ionic concentration x (here x = 0.5): (a) 
n = 0.7, U = 2, 4, and 6; (b) n = 0.3, U = 2, 3, 4, 5, and 6 
(U increases from bottom to top). Note the different range of 
A in both figures. 



We find a striking difference in the dependence of the 
Curie temperature T c on disorder strength A for differ- 
ent band fillings n < x and n > x (we chose x = 0.5 for 
numerical calculations). At n = 0.7, the critical tempera- 
ture T C (A) decreases with A for all values of U and even- 
tually vanishes at sufficiently large disorder [Fig. |l](a)]. 
By contrast, at n = 0.3, T C (A) weakly decreases with A 
at small U, but increases with A at large values of U 
[Fig. |(b)]. 

As will be explained below, this striking difference orig- 
inates from three distinct features of interacting electrons 
in the presence of binary alloy disorder: 

i) TP = T C (A = 0), the Curie temperature in the 
pure case, depends non-monotonically on band filling n. 
Namely, T£(n) has a maximum at some filling n = n*(U), 
which increases as U is increased jL7j; see Fig. |[ 

ii) In the alloy disordered system the band is split 
when A 3> W. As a consequence, for n < 2x and T« A 
electrons only occupy the lower alloy subband while the 
upper subband is empty. Effectively, one can therefore 
describe this system by a Hubbard model mapped onto 
the lower alloy subband. Hence, it corresponds to a single 
band with the effective filling n G ff = n/x. It is then 
possible to determine T c from the phase diagram of the 
Hubbard model without disorder [|l7| . 

iii) The disorder leads to a reduction of T^{n e q) by a 
factor x, i.e. we find 



T c (n) w xTP{n/x) 
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when A 3> W [^4|. Hence, as illustrated in Fig. ^, T c can 
be determined by T£(n e g). Surprisingly, then, it follows 
that, if U is sufficiently strong, the Curie temperature 
of a disordered system can be higher than that of the 
corresponding pure system [cf. Fig. ^|! 
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FIG. 2: Schematic plot explaining the filling dependence of 
T c for interacting electrons with strong binary alloy disor- 
der. Curves represent T£, the Curie temperature for the pure 
system, as a function of filling n at two different interactions 
Ui <C U'2 (cf. jTr[ ) . For n < x,T c of the disordered system can 
be obtained by transforming the open (for Ui) and the filled 
(for U2) point from n to n e ff, and then multiplying T£(n/x) 
by x as indicated by arrows. One finds T c (n) < Tf[n) for Ui, 
but T c (n) > Tc(n) for Ui- This difference originates from the 
non-monotonic dependence of T£ on n. 



To illustrate the alloy band splitting in the presence 
of strong interactions discussed above [see (ii)] we calcu- 
late the spectral density from the QMC results by the 
maximal entropy method |2l| . The results in Fig. ^ show 
the evolution of the spectral density in the paramagnetic 
phase at U = 4 and n — 0.3. At A = the lower and 
upper Hubbard subbands can be clearly identified. The 
quasiparticle resonance is merged with the lower Hub- 
bard subband due to the low filling of the band, and is 
reduced by the finite temperature. At A > the lower 
and upper alloy subbands begin to split off. A similar 
behavior was found at n = 0.7. The separation of the 
alloy subbands in the correlated electron system for in- 
creasing A is one of the preconditions [cf. (ii)] for the 
enhancement of T c by disorder when n < x, as discussed 
above. 

The splitting of the alloy subbands and, as a result, 
the changing of the band filling in the effective Hubbard 
model implies that T c vanishes for n > x. Namely, in the 
ferromagnetic ground state each of the alloy subbands 
can accommodate only xN a and (1 — x)N a electrons, re- 
spectively. Therefore, if the ground state of the system 
were ferromagnetic the upper alloy subband would be 
partially occupied for all n > x. This would, however, 
increase the energy of the system by A per particle in 
the upper alloy subband. Therefore, in the A 3> U limit 
the paramagnetic ground state is energetically favorable. 
This explains why T c vanishes at n = 0.7, as found in 
our QMC simulations [Fig. |l](a)]. Our conclusion that T c 
vanishes for n G ff = n/x > 1 when A 3> W is consistent 
with the observation in fun that there is no ferromag- 
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FIG. 3: Spectral density for different disorder strengths A 
at n — 0.3 and U = 4 as obtained by the maximal entropy 
method from QMC data at T = 0.071. The position of the 
lower/upper Hubbard subbands (LHB/UHB) are almost un- 
affected by the disorder, while the upper alloy subband shifts 
to the right as indicated by arrows. 



netism for n > 1 in the Hubbard model without disorder 
on fcc-lattice in infinite dimensions. 

The filling n — x is very particular because a new 
MIT of the Mott-Hubbard type occurs. Namely, when 
A increases (at U = 0), the non-interacting band splits, 
leaving 2xN a states in the lower and 2(1 — x)N a states 
in the upper alloy subband. Effectively, it means that 
at 11 = 1 the lower alloy subband is half filled (n c g = 
1). Consequently, a Mott-Hubbard MIT occurs in the 
lower alloy subband at sufficiently large interaction U 
p6f . In fact, for A ^> U we may infer a critical value 
U c = 1A7W* at T = from the results of Refs. [§7[ ||, 
where W* is the renormalized bandwidth of the lower 
alloy subband. Furthermore, from the analogy of this 
MIT with that in the pure case |29| we can expect a 
discontinuous transition for T < T* w 0.02VF*, and a 
smooth crossover for T > T* . From the results shown in 
Fig. I it follows that T* < 0.071, since for T = 0.071 and 
U = 6 a gap-like structure develops in the spectrum at 
A w 1.6, implying a smooth but rapid crossover from a 
metallic to an insulator-like phase p0| . 

The MIT described above is not obscured by the on- 
set of antiferromagnetic long-range order because in in- 
finite dimensions the fcc-lattice is completely frustrated 
p0| . Hence the insulator is paramagnetic. The transi- 
tion therefore occurs between a paramagnetic insulator 
(PI) at high T and a ferromagnetic metal (FM) at low 
T, at least at large U, as shown in the inset of Fig. |^. 
The actual boundary between the paramagnetic metal 
(PM) and the paramagnetic insulator-like phase has not 
yet been determined. The thick line in the inset of Fig. |] 
indicates the approximate position of the phase bound- 
ary between the PM and PI phases. We note that at the 
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FIG. 4: Spectral density for disorder strengths A = 0, 1, 1.6, 
and 1.8 (dotted, dashed, long-dashed and solid curves, re- 
spectively) at n = 0.5 and U — 6 as obtained by the maximal 
entropy method from QMC data at T — 0.071. For A > 1.6 a 
Mott-Hubbard gap-like structure develops around the Fermi 
level. Inset: A — T phase diagram of the binary alloy Hubbard 
model on the fcc-lattice in infinite dimensions at U — 6; PM - 
paramagnetic metal, PI - paramagnetic insulator-like phase, 
FM - ferromagnetic metal. Points with error bars represent 
the Curie temperatures obtained from QMC simulations; the 
solid line is a guide for the eye only. The thick line indicates 
the phase boundary between the PM and PI phases (see text). 
Circles: parameter values (A, T) corresponding to the spec- 
tral densities shown in the main panel. 



point where this boundary meets the FM phase the slope 
of the transition line T C (A) changes. 

In summary, we showed within DMFT that the inter- 
play between binary-alloy disorder and electronic corre- 
lation can result in unexpected effects, such as the en- 
hancement of the transition temperature T c for itiner- 
ant ferromagnetism by disorder, and the occurrence of a 
Mott-Hubbard type MIT off half-filling. An observation 
of these effects requires good control of the system pa- 
rameters over a wide range as was recently shown to be 
possible in experiments with optical lattices ]3l| . 
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